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Abstract.  Experiments  on  the  HMX-based,  condensed  explosive  PBX-9501  were  carried  out  to  validate  a 
reduced  asymptotically  derived  description  of  detonation  shock  dynamics  (DSD)  where  it  is  assumed  that 
the  normal  detonation  shock  speed  is  determined  by  the  total  shock  curvature.  The  passover  experiment 
has  an  embedded  lead  disk  in  a  right  circular  cylindrical  charge  of  PBX-9501  and  is  initiated  from  the 
bottom.  A  range  of  dynamically  changing  states,  with  both  divergent  (convex)  and  converging  (concave) 
shock  shapes  are  realized  as  the  detonation  passes  over  the  disk.  The  time  of  arrival  of  the  detonation  shock 
at  the  top  surface  of  the  charge  is  recorded  and  compared  against  the  DSD  simulation  and  a  separate  multi¬ 
material  simulation  (DNS).  A  new  wide-ranging  equation  of  state  (EOS)  and  rate  law  is  used  to  describe  the 
explosive  and  is  employed  in  both  the  theoretical  (DSD)  calculations  and  the  multi-material  simulations.  The 
experiment,  theory  and  simulation  are  found  to  be  in  excellent  agreement. 

Keywords:  Detonation  shock  dynamics,  PBX-9501,  experiment,  validation 


INTRODUCTION 

For  most  applications  of  condensed  explosives,  the 
length  of  the  reaction  zone  that  powers  the  lead  deto¬ 
nation  shock  is  a  fraction  of  a  millimeter  or  smaller. 
The  typical  ratio  of  the  device  size  to  the  reaction 
zone  is  typically  0(1000)  or  larger.  The  asymptotic 
theory  of  detonation  shock  dynamics  (DSD),  is  a  the¬ 
ory  that  assumes  that  the  radius  of  curvature  of  the 
shock  is  large  compared  to  the  length  of  the  reaction 
zone  that  supports  the  detonation  and  accounts  for 
change  to  the  detonation  shock  speed  due  to  shock 
curvature  during  quasi- steady  evolution.  The  present 
short  account  describes  a  validation  of  the  reduced 
theory  by  experiment  and  by  direct  simulation.  A 
full  account  of  this  work  and  complete  references  are 
found  in  [1],  and  [2]. 

The  simplest  DSD  theory  derives  the  result  that 
total  curvature  K  =  K\  +  K2,  is  a  function  of  the 
normal  detonation  shock  velocity,  Dn ,  written  as 


K  =  F(Dn ),  (1) 

with  the  property  F(Dcj)  =  0,  where  Dcj  is  the 
Chapman- Jouguet  velocity.  To  derive  (1)  from  the¬ 
ory,  one  assumes  that  the  explosive  is  described  by 
an  EOS  and  rate  law  of  the  general  form  e(p,v,X) 
and  r(p,  v,  A),  and  that  the  Euler  equations  hold.  One 
solves  for  the  quasi- steady  detonation  structure  in  the 
shock-attached  frame.  The  equations  for  the  struc¬ 
ture  contain  both  the  normal  detonation  speed  Dn 
and  the  total  curvature  K.  Since  it  is  assumed  that 
the  detonation  structure  passes  through  a  sonic  plane 
near  the  end  of  the  reaction  zone,  the  value  of  Dn  is 
not  independent  of  K  and  their  relationship  is  deter¬ 
mined  as  a  nonlinear  eigenvalue.  Since  the  normal 
shock,  curvature  relation  is  dependent  on  the  explo¬ 
sive’s  equation  of  state  and  reaction  rate  law,  experi¬ 
ments  serve  as  a  powerful  constraint  on  the  allowable 
forms  for  e(p,v,  A)  and  r(p,v,  A),  as  well  as  a  check 
on  the  theory. 
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FIGURE  1.  PBX  9501  Up  —  Us  Hugoniots  with  experi¬ 
mental  data  shown 

DESCRIPTION  OF  THE  EXPLOSIVE:  WIDE 
RANGING  EOS  AND  RATE  LAW 

Space  does  not  permit  recitation  of  the  detailed  forms 
used  for  the  equation  of  state  or  the  calibration  pro¬ 
cedure  used  to  assign  the  the  model  parameters  for 
the  wide  ranging  rate  law  for  PBX-9501;  full  details 
are  given  in  [1].  The  methods  used  in  the  calibra¬ 
tion  are  also  described  in  detail  in  [2].  Here  we  pro¬ 
vide  a  brief  summary.  Davis  developed  a  wide  rang¬ 
ing  equation  of  state  for  detonation  products  whose 
form  was  chosen  to  accurately  describe  the  behav¬ 
ior  of  adiabatic  y  (dimensionless  sound  speed)  and 
Griineisen  gamma,  F.  Davis  also  developed  a  similar 
reactants  equations  of  state  fit  to  PBX-9404,  9501. 
Stewart,  Yoo  and  Davis  in  the  12-th  Detonation  Sym¬ 
posium,  proposed  a  modification  of  Davis’  reactant 
equation  of  state  and  introduced  a  closure  model  to 
develop  a  mixture  EOS  that  includes  the  reaction 
progress  variable  of  the  form  e(p,v.  A)  and  that  uses 
the  standard  rules  for  a  binary  mixture  of  reactants 
and  products,  where  A  is  the  mass  fraction  of  the 
products  and  is  the  reaction  progress  variable.  The 
equation  of  state  parameters  were  fit  to  the  shock 
Hugoniot  data  for  both  reactants  and  products.  Fig. 
1  shows  a  plot  of  the  particle  velocity,  shock  veloc¬ 
ity  (Up,  Us)  Hugoniot  (top  curve)  calculated  from  the 
products  EOS,  and  Hugoniot  (bottom  curve)  calcu¬ 
lated  from  the  reactants  EOS,  compared  with  experi¬ 
ment.  The  experiment  data  shown  was  was  compiled 
by  R.Gustavsen,  LANL.  The  calibration  considers 
work  done  by  expanding  gases  and  the  temperatures 
of  the  reactants  and  products. 

A  rate  law  is  proposed  for  PBX-9501  with  a 
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FIGURE  2.  "Pop "-plot,  run  to  detonation  distance  ver¬ 
sus  input  shock  pressure  (Pinput)  for  PBX  9501  from  exper¬ 
iments  and  direct  simulation  results. 

single-term,  fractional  depletion,  pressure  dependent 
reaction  rate  of  the  form 

r(p,v,h)  =k(\-X)v  .  (2) 

Shock  to  detonation  data  ("Pop" -plot)  and  detona¬ 
tion  shock  speed  curvature  data,  is  used  to  calibrate 
the  parameters  of  the  rate  law.  Hull’s  experimental 
data  suggests  the  Dn ,  K  relation  is  linear  near  the  CJ 
point.  The  depletion  exponent  v  is  picked  primarily 
to  match  Hull’s  data,  see  Fig.  3.  The  pressure  expo¬ 
nent  N  and  rate  constant  k  are  adjusted  to  match  the 
shock  initiation  data.  One-dimensional,  reverse  im¬ 
pact  simulations  were  carried  out  using  the  specified 
EOS  and  rate  law  to  match  the  published  experimen¬ 
tal  data  for  PBX-9501,  and  the  results  are  shown  in 
Fig.  2.  The  calibrated  rate  law  parameters  were  found 
to  be  pcj  =  36.3  GPa,  k  =  110 psec~l,  N  =  3.5  , 
v  =  0.93. 

Figure  3  shows  the  detonation  velocity  curvature 
relation  calculated  using  this  EOS/rate  law  pair  de¬ 
scribed  by  asymptotic  theory.  This  Dn  —  K  relation 
is  the  shock  motion  rule  that  is  used  to  compute  the 
shock  motion  according  to  the  reduced  DSD  descrip¬ 
tion.  Hull’s  Dn ,  K  experimental  data  is  also  shown. 


THE  PASSOVER  EXPERIMENT 

The  experimental  set  up  is  shown  in  Fig.  4.  The 
PBX-9501  explosive  (white  material)  has  a  disk  of 
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FIGURE  3.  The  Dn ,  K  relation  for  PBX-9501  calculated 
from  the  wide-ranging  EOS  and  rate  law  model 


pure  lead  (grey  object)  embedded  along  the  cen¬ 
tral  axis.  A  detailed  description  of  the  experiment 
is  found  in  [1].  The  experiment  transforms  a  sin¬ 
gle,  quasi-steady,  convex  hemispherical  shock  into 
a  shock  with  a  high  concavity  at  the  central  implo¬ 
sion  axis.  PBX-9501  was  chosen  to  test  previous  ex¬ 
perimental  DSD-characterization  from  Hull’s  work. 
Pure  lead  was  selected  as  the  inert  for  its  high  shock- 
impedance  and  well-characterized  shock  Hugoniot 
properties. The  charge  is  initiated  at  the  bottom  of 
the  PBX-9501  charge.  The  detonation  shock  front 
propagates  as  a  simply  connected  surface  with  con¬ 
vex,  positive  curvature  hemisphere  for  50-mm  and 
then  encounters  the  lead  disk  within  the  top  piece 
of  PBX-9501.  The  shock  speed  in  the  inert  lead  is 
much  lower  than  the  detonation  velocity  in  the  PBX- 
9501  and  a  diffraction  event  occurs  as  the  detona¬ 
tion  sweeps  about  the  disk  and  encompasses  it.  Wa¬ 
ter  sits  atop  the  charge  and  extinguishes  the  reactive 
shock  as  it  transmits  into  the  water.  Light  from  the 
explosive/water  plane  records  the  time-of-arrival  of 
the  detonation  shock  through  a  single  150-micron  slit 
aperture  plate  as  captured  by  a  Cordin  132A  cam¬ 
era.  Four  (4)  passover  experiments  were  conducted 
with  identical  hardware.  Figure  5  (as  computed  from 
DSD-theory)  illustrates  the  subsequent  shock  motion 
in  the  experiment. 

COMPARISON  OF  DSD,  NUMERICAL 
SIMULATIONS  (DNS)  AND  EXPERIMENT 

Two  different  types  of  simulations  were  carried  out 
to  compare  with  the  time  of  arrival  results  obtained 


FIGURE  4.  Assembly  sketch  for  DSD  validation  exper¬ 
iment 


FIGURE  5.  DSD  simulation  of  the  axisymmetric 
passover  experiment.  The  grey-scale  shows  the  shock  pres¬ 
sure  in  GPa  when  a  shock  passes  a  point  (x,y)  in  the  ex¬ 
plosive. 


by  the  passover  experiments.  The  first  uses  the  re¬ 
duced  DSD-model  defined  by  the  Dn ,  K  shock  mo¬ 
tion  rule,  shown  in  Fig.  1,  subject  to  inert  angle  con¬ 
finement  boundary  conditions.  The  lead  disk  is  de¬ 
scribed  by  a  Mie-Gruneisen  (UP,US)  EOS  e(p,v)  . 
Shock  polar  analysis  for  the  DSD  confinement  de¬ 
rives  the  angles  at  the  PBX-950 1/lead  interface  (the 
interior  angle  between  the  shock  and  interface  nor¬ 
mals)  as  co  =  35° (sonic)  and  66° (subsonic).  The  ini¬ 
tial  shock  in  the  DSD  simulation  was  a  hemi-sphere 
of  radius  5  mm  centered  at  the  bottom  of  the  charge. 

The  multi-material  numerical  simulation  (DNS) 
was  carried  out  with  the  wide-ranging  EOS  and  rate 
law  for  the  explosive  and  the  Mie-Gruneisen  EOS  for 


the  lead.  An  outflow  boundary  condition  is  used  at 
the  lateral  boundaries.  The  multi-material  simulation 
code  combines  two  high-order  solvers,  a  high-order 
total  variation  diminishing  (TVD)  solver  for  the  Eu¬ 
ler  equations  and  a  level-set  solver  to  move  the  ma¬ 
terial  interface  that  separates  explosive  and  lead.  The 
DNS  simulation  uses  an  initial  condition  of  a  hemi¬ 
spherical  hot  spot  of  radius  R  =  5  mm,  centered  at  the 
bottom.  Figure  6  shows  a  comparison  of  the  DSD  and 
DNS  simulation,  to  show  that  they  give  consistent  re¬ 
sults  (i.e.  that  the  shocks  overlap). 


CONCLUSION 

Figure  7  is  a  composite  of  TOA  records  for  the 
passover  experiments  that  includes  experiments, 
DSD  simulations,  DNS  simulations  and  the  ideal 
Huygens  construction.  The  Huygens  construction 
assumes  that  the  detonation  propagates  at  a  constant 
normal  shock  speed  of  Dcj  =  8.86  mm/flsec  and 
does  not  account  for  the  slowing  of  the  wave  due  to 
curvature  effects.  The  excellent  level  of  agreement, 
both  qualitative  and  quantitative,  between  exper¬ 
iment,  DSD  and  DNS  is  encouraging  because  it 
indicates  that  one  can  use  the  wide  ranging  EOS/rate 
law  and  the  corresponding  DSD  description  effec¬ 
tively  to  model  real  explosives  and  predict  complex 
dynamic  behaviors. 
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FIGURE  6.  Comparison  of  DSD  and  DNS  simulation 
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FIGURE  7.  TOA  comparison  of  experiments,  DNS,  and 
DSD  simulations  at  the  top  of  the  charge 


